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Abstract 


Surface albedo is a key parameter in the Earth’s energy balance since it affects the amount of solar 
radiation directly absorbed at the planet surface. Its variability in time and space can be globally 
retrieved through the use of remote sensing products. To evaluate and improve the quality of satellite 
retrievals, careful intercomparisons with in situ measurements of surface albedo are crucial. For this 
purpose we compared MODIS albedo retrievals with surface measurements taken at 53 FLUXNET 
sites that met strict conditions of land cover homogeneity. A good agreement between mean yearly 
values of satellite retrievals and in situ measurements was found (R 2 = 0.82). The mismatch is 
correlated to the spatial heterogeneity of surface albedo, stressing the relevance of land cover 
homogeneity when comparing point to pixel data. When the seasonal patterns of MODIS albedo is 
considered for different plant functional types, the match with surface observation is extremely good at 
all forest sites. On the contrary, in non-forest sites satellite retrievals underestimate in situ 
measurements across the seasonal cycle. The mismatch observed at grasslands and croplands sites is 
likely due to the extreme fragmentation of these landscapes , as confirmed by geostatistical attributes 
derived from high resolution scenes. 

Keyword: MODIS, surface albedo, validation, FFUXNET, terrestrial ecosystems, plant functional 
types, remote sensing, 

1. Introduction 

Fand surface broadband albedo directly affects Earth’s climate by determining the fraction of 
shortwave radiation absorbed at the ground and therefore influencing the surface energy budget 
(Dickinson, 1983). Surface albedo is a crucial parameter in determining the magnitude of energy fluxes 
in the soil-plant-atmosphere continuum (Bonan, 2008; Chapin et al., 2008), affecting surface 


temperature, evaporation and transpiration, cloud formation and precipitation, thus ultimately 
impacting gross primary productivity (Dickinson, 1983; Fawrence & Slingo, 2004; Ollinger et al., 
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2008; Sellers et al., 1997). Several authors have investigated the interplay between albedo and drought 
(Govaerts & Lattanzio, 2008) or fires (Randerson et al., 2006), and the climate sensitivity to variation 
in surface albedo caused by major changes in land cover as the expansion of agricultural land in the 
northern hemisphere during the 18 th century (Myhre et al., 2005; Vavrus et al., 2008). Surface albedo is 
also a key factor in the expected positive feedback between surface temperature and global warming at 
northern latitudes (Chapin et al., 2005) and may play a relevant role in offsetting the carbon 
sequestration potential of afforestation programs (Anderson et al., 2010; Betts, 2000; Betts et al., 2007; 
Bird et al., 2008; Rotenberg & Yakir, 2010). 

Given the relevance of surface albedo in the Earth’s climate system, monitoring this parameter in space 
and time is fundamental for the development of global climate models (Alton, 2009; Frida A-M et al., 
2006; Hollinger et al., 2009; Tian et al., 2004) and for climate change and ecosystem research in 
general (Betts, 2000; Charlson et al., 2007; Chamey et al., 1977; Dirmeyer & Shukla, 1994; Hall & Qu, 
2006; Henderson-Sellers & Wilson, 1983; Pinty et al., 2011a). An important step toward the 
availability of global surface spectral albedo has been the launch of NASA’s Terra and Aqua satellites 
and the MODerate-resolution Imaging Spectro radiometer (MODIS) (Lucht et al., 2000b; Salomonson 
et al., 1989; Schaaf et al., 2002). The MODIS sensor provides global maps of surface albedo 
reconstructed from retrieved models of reflectance anisotropy at a 500-m gridded spatial resolution 
every 16 days for the first seven MODIS spectral bands (0.47-2. 1 pm) and for three broadband regions 
(0.3-0. 7, 0. 7-5.0, and 0. 3-5.0, pm) (Lucht et al., 2000b; Moody et al., 2008; Schaaf et al., 2002). 
Comparing satellite albedo retrievals with surface measurements and with independent satellite 
products is fundamental in evaluating the accuracy of remote sensing products and improving retrieval 
algorithms (Liang et al., 2002; Pinty et al., 2011b). Several recent studies have evaluated the 
consistency of global albedo products using in situ data at various spatial and temporal scales (Chen et 
al., 2008; Jin et al., 2003a; Jin et al., 2003b; Liang et al., 2002; Liu et al., 2009; Roman et al., 2010; 

Roman et al., 2009; Wang et al., 2010) and under specific snow cover conditions (Stroeve et al., 2005). 
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Most of these studies stress that a direct comparison is very challenging because of scale mismatch and 
heterogeneity of the land surface at the satellite measurement scale that reduces the spatial 
representativeness of ground point measurements (Liang et al., 2002; Roman et al., 2010; Roman et al., 
2009). As a consequence, a careful selection of ground points and the characterisation of their spatial 
representativeness are crucial for a meaningful point-to-pixel comparison (Liang et al., 2002; Lucht et 
al., 2000a; Roman et al., 2009). 

Intercomparisons of surface and satellite albedo have been performed so far at a limited number of 
locations (Jin et al., 2003a; Liu et al., 2009; Roman et al., 2010; Roman et al., 2009; Salomon et al., 
2006; Wang et al., 2010) and a global analysis across different continents and plant functional types 
(PFTs) is still lacking. The objective of this work is to provide a comprehensive intercomparison in 
time and space of in situ measurements and satellite retrievals of snow-free broadband surface albedo. 
For this purpose we compared MODIS gridded albedo retrievals at the 500-m scale with ground 
measurements performed across the FLUXNET network (Baldocchi et al., 2001), the largest global 
data set of energy and mass flux measurements at ecosystem scale. 

The geographical extent of the terrestrial data set allowed the comparison of several PFTs in a 
comprehensive and consistent way across the seasonal cycle. In addition, the large number of 
experimental sites in the network provided an unprecedented opportunity to perform a careful 
evaluation of the surface heterogeneity at the reference plots, based on a combination of qualitative and 
quantitative metrics. For this purpose images from Google Earth, MODIS and Enhanced Thematic 
Mapper Plus (ETM+) have been used at various spatial scales (from lxl to 7x7 km). Differences 
between satellite retrievals and in-situ albedo have been analyzed as a function of surface 
heterogeneity, PFT and seasonality. Results of the intercomparison have been finally discussed 
considering the different sources of uncertainty that affect the terrestrial and satellite datasets. 
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2. Materials and methods 


2.1. Surface data set 

In this study, we used in situ radiometric measurements available in the FLUXNET “La Thuile” 
database (www.fluxdata.org , October 2010) released in December 2007, which includes half hourly 
observations of ecosystem fluxes and meteorological data from more than 250 sites, for a total of 960 
site-years. 

Albedo is computed as the ratio of downward and upward global radiation as observed with double 
pyranometers (e.g. CMA-1 1, CMA-6 or CNR-1, Kipp&Zonen, Delft, The Netherlands). Surface albedo 
is typically estimated in the spectral range 280-2800 nm (accounting for more than 98.5% of the 
surface solar radiation according to ASTM G-173 reference spectra) and is therefore comparable with 
the broadband MODIS albedo (300-5000nm). Giving that the field of view (FOV) of pyranometers is 
typically 180 deg, the footprint of surface reflectance measurements is theoretical infinite. Flowever, 
due to the cosine response of the sensor, 50% of the signal originates in a FOV of 90 deg and 80% in a 
FOV of 127 deg. The footprint of surface albedo therefore depends on the height of the albedometer 
above the canopy top (ranging from 5 to 10 m) and typically extends up to 10-20 m from the tower at 
80% of the signal. 

The uncertainty of surface albedo measurements depends on the absolute accuracy of phyranometers 
(about 5%) and on the non-ideal cosine response (about 3%). Most of the errors associated with the 
absolute accuracy of the instrument are similar for upward and downward fluxes and therefore 
compensate. Overall the expected accuracy is in the order of 4-7% in clear sky and 1-4% in overcast 
condition (Pirazzini, 2004; Pirazzini et al., 2006). 

The geographical distribution of the sites is strongly clustered in Europe and North America (97 and 
106 sites corresponding to 38% and 42% of the total), which are the regions with the longest history of 
continuous ecosystem flux measurements (Baldocchi et al., 2001). Several sites in the database are 
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located in tropical Amazonia and East Asia, while the coverage in Africa, Central Asia, and Australia 
remains sparse and limited in the number of observation years. Despite the uneven geographical 
distribution, the “La Thuile” database guarantees a good coverage of the most important plant 
functional types, among which evergreen needle leaf forest (ENF), grassland (GRA), deciduous 
broadleaf forest (DBF), and cropland (CRO) are the most represented with respectively 28%, 18%, 

13% and 12% of the sites. 

Out of the 138 FLUXNET sites reporting continuous measurements of incoming and outgoing 
shortwave radiation (300-2800 nm; CMA-1 1, CMA-6 or CNR-1, Kipp&Zonen, Delft, The 
Netherlands) 18 have been excluded after a QA/QC analysis of the albedo data series. The QA/QC 
procedure was based on the following criteria: occurrence of an offset in the incoming or reflected 
radiation (night-time data systematically and significantly different from zero), occurrence of phase lag 
between incident and reflected radiation and systematic occurrence of unrealistic values (e.g. reflected 
radiation higher than incident radiation). 

The land cover characteristics of the remaining 120 sites have been carefully classified using high 
resolution satellite images (available via Google Earth™ ), to identify those matching the requirement 
of homogeneity in the area surrounding the measurement tower (Jin et al., 2003b; Roman et al., 2010; 
Roman et al., 2009). Although MODIS albedo is gridded at 500-m resolution the land classification has 
been performed at 1km 2 , taking into account the uncertainty in the geospatial registration of satellite 
products and the fact that the albedo retrieval algorithm is based on multi-angle observations covering 
larger areas at edge of scan. 

The classification process was based on the following four steps: 

1. visual identification of the number and extension of different PFTs in the 1km 2 area 
surrounding the tower; 

2. verification of the correspondence between the dominant PFT in the 1km 2 area and the PFT at 
the tower site as reported in the FLUXNET database; 
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3. qualitative ranking of landscape heterogeneity in three classes (low, medium, high) based on the 
plant canopy characteristics (tree density, patchiness, etc.); 

4. attribution of a confidence level in the classification of the sites (low, medium, high) based on 
the quality of the image. 

To guarantee the highest level of homogeneity and to minimize issues associated with spatial 
representativeness in the point-to-pixel comparison, only those sites characterized by the lowest level 
of heterogeneity and with only one PFT in the 1km 2 area were included in the analysis. 

2.2. MODIS products 

The MODIS albedo retrievals at the FLUXNET sites were generated using three MODIS products, 
namely, MCD43A1 (BRDF-Albedo Model Parameters 16-Day L3 Global 500m), M*D04 (Aerosol 
product daily L2 Global 10km), and MCD43A2 (BRDF-Albedo Quality 16-Day L3 Global 500m). All 
these products are from the Collection V005 MODIS reprocessing campaign. The MODIS surface 
reflectance anisotropy and albedo product is based on all high quality, cloud-free, atmospherically 
corrected surface reflectances that are obtained over a 16-day period. When sufficient observations are 
available to adequately sample the surface anisotropy, an appropriate rendition of the 
RossThickLiSparseReciprocal Bidirectional Distribution Reflectance Model (BRDF) model is retrieved 
(Lucht et al., 2000b; Schaaf et al., 2002). This retrieval is attempted every 8 days at a 500 m gridded 
resolution. This retrieval model is used to generate intrinsic values of clear-sky direct surface albedo 
(referred to as directional hemispherical reflectance or black-sky albedo) and wholly diffuse albedo 
under isotropic illumination (bihemispherical albedo or white-sky albedo). These can be combined 
under particular illumination and atmospheric aerosol optical depth conditions (Lucht et al., 2000b; 
Roman et al., 2010) to provide clear-sky albedos comparable to those measured in situ at a flux tower. 
Albedo quantities are reported at a 500-m gridded resolution, but all multi-angle observations that 
encompass areas are utilized in the retrieval. Therefore, although extended observation coverage is 
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somewhat compensated for in the retrieval process, it is best to consider regions larger than 500 m 
when comparing observations made from satellite to those made on the ground. 

The calculation of clear-sky surface albedo at the tower sites involved the following two steps. The first 
step was the generation of the aerosol optical depth values for each site and each calendar date using 
the MODIS - Terra (MOD04) and MODIS-Aqua (MYD04) aerosol swath products. To generate the 
optical depth, the MODIS Adaptive Processing System - MODAPS (Masuoka et al., 2000; Masuoka et 
al., 2007) was used to prepare M*04 subsets at 50 x 50 km region centered at the site. All pixels that 
had optical depth values greater than 0.35 or a cloud fraction greater than 0.6 were filtered out and not 
used in the optical depth generation. All pixels that had fill values for solar zenith angle were also 
rejected. After the filters were applied, a combined M[OY]D optical depth file was generated for each 
site, taking valid optical depth values from Terra and Aqua and generating one mean value for the 
optical depth per site per day. This method of course is not as accurate as having instantaneous sun 
photometer data (Holben et al., 1998) at the site, but the mean gives an approximation of the aerosol 
optical depth over the local solar noon. 

The second step was the calculation of the clear sky surface albedo on the basis of the MODIS-derived 
550nm aerosol optical depths calculated in the previous step, the local solar zenith angle, the 
MCD43A1 product, and QA flags from MCD43A2 for each site involved in the analysis and for each 
date. If a date has no valid MCD43A1 pixels or if the optical depth was a fill value, no albedo was 
calculated for that date. As far as quality criteria are concerned only “full BRDF inversion” pixels 
(QA=0 processed, good quality) were included in the calculation, while the "Snow BRDF Albedo" 
band of the MCD43A2 product was used to identify and exclude snow albedo retrievals. Following this 
procedure clear-sky MODIS albedo at local solar noon were retrieved at each FLUXNET site for all 
days with available aerosol MODIS product (M*D04) information, snow-free conditions, and solar 
elevation angles greater than 20 deg. On the same dates, the flux tower measurements of albedo have 
been averaged for the hour centered at solar noon. 
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To integrate the observations at the FLUXNET sites in the global picture, snow-free global albedo 


averages per PFT and latitudinal band were computed from the MODIS V005 0.05 degree Climate 
Modeling Grid (CMG) product and stratified with the MCD12C1 land cover product. Yearly averages 
have been calculated on each pixel fulfilling the following requirements: QA=0 (majority processed, 
good quality), snow coverage less than 10% (based on MODIS estimates), and major PFT coverage 
greater than 70% of the pixel. Note that the 0.05 degree MCD43C1 product is an average of the 500m 
pixel underlying each 0.05 degree pixel and the quality flag only represents the quality of the majority 
of the underlying pixels. 

2.3. Landscape heterogeneity 

One of the key issues in the intercomparison of satellite retrieval and surface observations is the 
objective and quantitative evaluation of landscape heterogeneity and the representativeness of in situ 
measurements (Liang et al., 2002; Roman et ah, 2009; Susaki et ah, 2007). 

For this purpose we applied the methodology presented by Roman et al. (2009) and based on the 
estimation of geostatistical attributes from high resolution scenes (Enhanced Thematic Mapper Plus). 
The spatial patterns and scales of landscape heterogeneity have been estimated from variogram models 
fitted at FLUXNET sites over the spatial scales of MODIS observations. 

In synthesis, the methodology adopted for the estimation of geostatistical indexes is based on the 
comparison of variogram model parameters retrieved at different spatial resolution (i.e. from 1.0 km 2 to 
1.5 km 2 squared subsets). By examining the variogram parameters at two scales, the spatial 
characteristics of a given measurement site is compared against the larger landscapes extending to 
several MODIS pixels. 

Four different geostatistical attributes of spatial representativeness have been used to describe the 
overall variability (Rev), spatial extent (Rse), strength of the spatial correlation (Rst), and 
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1 spatial structure (Rst) of surface albedo for a given measurement site. Further details on the methods 

2 and algorithms used to calculate these attributes are reported in Roman et al. (2009). 

3 Using a weighted combination of the four geostatistical attributes a comprehensive metric of the 

4 landscape heterogeneity (ST scor e) has been computed to evaluate and compare FLUXNET sites: 

IR I +IR l + IR I 

score ^ 

5 In addition, the landscape heterogeneity at a larger spatial scale (7x7 km) has been quantified as the 

6 standard deviation of MODIS albedo (14 x 14, 196 pixels) area centered at the FLUXNET site. 
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3. Results and discussion 


3.1. The terrestrial data set: spatio-temporal distribution and representativeness 

As a result of the visual classification of land cover characteristics, 49% of the 120 sites performing 
reliable continuous measurements of broadband albedo were rejected for reasons of landscape 
heterogeneity and 7% because of low confidence in the PFT classification. 

Of the remaining 53 sites (among which there were 15 ENF, 8 CRO, 7 DBF, and 7 evergreen broad- 
leaf forests (EBF)) the largest fraction is located in Europe (19) and North America (15 USA, 7 
Canada). The remainder of the sites are located in Africa (4), South America (3), Asia (3) and Australia 
(2) (Fig. 1, Table 1). 

[Figure 1] 

[Table 1] 

Across all sites, 18666 days of synchronous MODIS retrievals and in situ surface albedo were 
available. On average, this worked out to 333 and 400 days of data for each forested and non-forested 
site, respectively (on average about 80 observations per site and year), but there was considerable 
variability among sites. To date, this is the largest data set that has been used to compare satellite and in 
situ albedo measurements. 

Although the data set spans the years 2000-2007, the vast majority (87%) of observations were made 
between 2003 and 2006 (Fig. 2a). Because of filtering for snow and cloud cover, and the dominance of 
Northern hemisphere sites in the data set, the number of observations is lower during winter months 
(November through February) (Fig. 2b). 

[Figure 2] 

The seasonal variation of the global average of snow-free albedo is remarkably small for forests, likely 
because ~70% of the sites are evergreen forests with leaf area index, canopy structure and chemistry 


not as dynamic as those of deciduous ecosystems. The seasonal trend of snow-free albedo is somewhat 
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more pronounced for non-forest ecosystems, which include crops and grasslands (Running et al, 1995) 
(Fig. 2b). 

The spatial representativeness of the FLUXNET albedo data set in the global albedo domain has been 
explored by superimposing the MODIS retrievals at the sites on the global distribution of MODIS 
albedo (Fig. 3) for the same PFT or latitudinal class. Concerning the latitudinal distribution, 75% of 
FLUXNET sites are clustered in the 30-55° N band (Fig. 3a), and in this latitude band the surface 
albedo observed at the sites is substantially lower (20%) than the MODIS latitudinal average. This is 
due to the non-representative distribution of FLUXNET sites, which are mostly located in dense and 
productive temperate forests of the Northern Flemisphere. The difference between latitudinal MODIS 
averages and site-measurements of albedo is less evident in the Southern hemisphere, where also the 
number of FLUXNET site is remarkably lower. 

[Figure 3] 

In terms of PFTs (Fig 3b), 79% of sites analyzed here are classified as one of ENF, EBF, CRO, GRA or 
DBF. When PFT averages are compared with the global average of MODIS albedo, most of the sites 
fall within the 10-90 percentile intervals of the global MODIS observations, with the exception of EBF 
for which MODIS retrievals at the sites show lower values of albedo. For the other PFTs (woody 
savanna (WSA), open shrubland (OSF1), savanna (SAV), mixed forest (MF)) the number of sites is too 
low to speculate on the global representativeness of the FLUXNET dataset. 

3.2. Point to pixel comparison 

The comparison of MODIS albedo retrievals with in situ measurements has been limited to the 
FLUXNET sites with the highest degree of homogeneity in order to minimize the effect of the scale 
mismatch. Examples of high resolution images (available via Google Earth™) used to evaluate site 
homogeneity are reported in Fig. 4 for four test sites. 
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[Figure 4 ] 

The time series of MODIS retrievals and in situ measurements at these four sites show contrasting 
results (Fig. 5). At some sites the match is extremely good both in terms of absolute values and 
seasonal trend (e.g. US-MMS, US-FPe). At other locations MODIS retrievals show a systematic 
overestimation (e.g. PT-Esp) or underestimation (e.g. ES-LMa) of surface measurements. These biases 
are probably due to the fine-scale spatial variability of the plant cover and to the representativeness of 
the tower footprint in the MODIS pixel, since the other sources of uncertainties (sensor calibration, 
uncertainty of AOD estimates, etc.) cannot explain such large and systematic errors. In particular at 
sites with discontinuous plant canopies like ES-LMa (Fig. 4), the height from the ground and the spatial 
location of the albedometer are critical factors, determining the representativeness of the in situ 
measurements at the resolution of the satellite pixel. 

[Figure 5] 

A good agreement between satellite and in situ measurements was found when sites were grouped 
according to PFT (Fig. 6a, type II regression, R 2 = 0.82). Forest PFTs (such as ENF, EBF, MF) fall in 
the lower part of the graph, while non-forest PFTs (such as CRO, GRA, wetlands (WET)) in the upper 
part. Albedo is correlated to several plant-level traits including leaf albedo, leaf area index, vertical 
angle of leaves/needles, the degree of foliage and canopy clumping, and the geometric-optical 
shadowing due to canopy structure. It is well known that the structural canopy traits typical of tall 
canopies trap more of the incoming radiation, therefore reducing the canopy albedo of forest PFTs 
(Cescatti, 1998; Davidson & Wang, 2004). 

[Figure 6] 

Satellite retrievals and surface measurements at the site level (Fig. 6b) show a coefficient of 
determination of 0.83, very similar to that observed in the comparison by PFT and reported in 
comparable studies (Wang et al. 2010). Both regressions in Fig. 6 show a slope lower than one (0.88 
and 0.74 for panel (a) and (b), respectively), primarily due to the different spatial scale of the two 
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estimates, with the satellite retrievals referring to a considerably larger area than tower measurements. 
The spatial averaging of surface albedo on the large MODIS pixels dampen the variability when 
compared with high-resolution imagery or point measurements and ultimately determines the trend 
(intercept >0 and slope <1) of the regressions in Fig. 6. Similar slope and closer correlations have been 
observed by comparing early MODIS retrievals with surface measurements using high-resolution 
remotely sensed imagery (Landsat7 ETM+) to characterize the land cover heterogeneity in the MODIS 
pixel (Liang et al., 2002). The negative bias of MODIS at the sites with larger albedo has been reported 
also by Wang et al. (2010) and has been mostly attributed to the underestimation of visible albedo. 

The mean error and mean absolute error of the MODIS retrievals vs. the in situ observations are 
reported in Fig. 7a together with the R 2 of the regression computed for the single sites. The mean error 
is a measure of the retrieval accuracy and is on average negative (-0.004) implying that the magnitude 
of the MODIS surface albedo is on average slightly lower than in situ observations. It is also interesting 
to notice that the spread is considerably larger among the sites with a negative error both in terms of 
mean error and mean absolute error. The R 2 observed at single sites, proportional to dot size in Fig. 7a, 
is largely independent of the mean error since it mostly depends on the seasonal variability of the 
measurements. For this reason, sites with a marked seasonality in albedo (i.e. deciduous forests, crops 
and boreal ecosystems) typically show a higher R 2 than tropical or Mediterranean evergreen forests. 
The error distribution peaks in the -0.01 class and is skewed to the left, with larger errors at sites where 
MODIS underestimates surface measurements (Fig. 7b). The mean values of STscore (geostatistical 
index proportional to the representativity of in-situ measurements in the MODIS pixel) in the different 
classes of mean error shows several interesting features (Fig. 7b). The highest mean value of the score 
(3.9) is in the class at zero mean error, confirming the validity of this approach to quantify the spatial 
representativeness of experimental sites. The values of ST scor e is considerably lower (1-1.4) at sites 
characterized by a negative error (MODIS<in situ) and the variability of this index is rather limited 
between sites. At these sites the towers are probably located in a brighter spot than the average of the 


15 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 


pixel and typically pertain to the categories of grassland (GRA) or cropland (CRO) (Table 2). On the 
contrary sites with a positive bias (MODIS>In situ albedo), where the tower spot is located on an area 
darker than the surroundings, show intermediate value of the score (2.5) and a large variability between 
sites. 

[Figure 7] 

[Table 2] 

Giving that the MODIS albedo product is reported at 500 m resolution while the retrieval algorithm is 
based on multiangular reflectances that extend over this area, landscape heterogeneity at scales larger 
than 500 m may ultimately affect the uncertainty in the retrievals. To assess the impact of the spatial 
heterogeneity on the accuracy of MODIS retrievals we computed the mean absolute percentage error 
for sites ensembles at increasing levels of spatial variability (expressed as standard deviation of 
MODIS albedo in the 7x7 km 2 area surrounding the site or as ST SCO re). Figure 8 shows that the error 
between satellite and in situ measurements increases with the spatial variability of albedo, further 
demonstrating the crucial issue of the spatial variability of surface albedo in the point to pixel 
intercomparisons. 

Both MODIS standard deviation and ST score show sharp variations for values of the mean absolute 
percentage error up to 12%. Above this threshold both indexes are no longer correlated with the error. 
The trend of the two spatial indexes is remarkably similar though in opposition, since MODIS standard 
deviation increases with the spatial variability while ST score increases with spatial homogeneity. Despite 
the different spatial domain (7x7 and 1.5 km) the two indexes are strongly correlated (Fig. 8 inset) with 
an R 2 of 0.90. These results demonstrate that a close match between surface measurements and satellite 
retrievals is achievable only at the most spatially homogenous sites, where the canopy optical 
properties are scale-invariant. These golden sites can be effectively detected with both statistical 
analysis of MODIS albedo or with higher resolution images and detailed geostatistical indexes. 
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[Figure 8] 


3.3. Seasonal trends 

Seasonal patterns in the relationship between site-level averages of MODIS and FLUXNET albedo are 
reported in Figure 9. In all seasons, the relationship between MODIS (y axis) and in situ (x axis) albedo 
has a slope ranging from 0.71 to 0.74 and positive intercept, meaning that MODIS over-estimates 
surface albedo of low-albedo sites, and under-estimates albedo of high-albedo sites. The lowest 
correlation (R 2 =0.73) is observed for winter values (from December to February) and the highest 
(R =0.82) for fall values (from September to November). Similarly, larger differences between in situ 
and satellite based estimation of surface albedo in winter season were observed in Jin et al. (2003), who 
suggested that this was the result of the increased heterogeneity of surface reflectivity due to the 
presence of residual snow and canopy heterogeneity. 

[Figure 9] 

The inset in Fig. 9a shows the correlation matrix of the errors observed in the different seasons. All 
correlation coefficients are positive, meaning that sites with surface measurements larger than satellite 
retrievals in one season tend to have the same behaviour in the other seasons. Thus, site-specific biases 
are coherent over time, possibly indicating mismatches between the footprint of tower radiometric 
instruments and the corresponding MODIS pixel, calibration errors of the tower radiometric 
instruments, inappropriate estimation of aerosol optical depths or systematic errors in MODIS retrieval. 
[Figure 10] 

In Figure 10 the seasonal trends of satellite and in situ data for a set of PFTs are reported with their 
standard deviation, and superimposed on the range of MODIS albedo observed at the global scale for 
the same PFTs (grey bands represent the 10-90 percentile interval). An exceptionally good agreement 
between the two independent albedo estimates is observed for forest PFTs, both in terms of absolute 
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albedo values and seasonal patterns. In contrast, for non-forest PFTs (CRO and GRA) MODIS 
retrievals are systematically smaller than surface measurements, while their respective seasonal 
patterns are in agreement as observed also by Davidson & Wang (2004). The mismatch observed at the 
non-forest sites is likely due to the extreme fragmentation of these landscapes and to the very limited 
spatial footprint of surface radiometric measurements performed at these sites (i.e. over short 
vegetation albedometers are typically installed few meters above ground). Wang et al. (2010) reported 
that the average negative bias of MODIS albedo is probably due to the underestimation of visible 
surface reflectance and that this phenomenon may be more pronounced for herbaceous canopies with 
larger albedo in the visible bands. 

The large variability of CRO and GRA surface albedo (Fig. 10) can be ascribed to the spatial variability 
in LAI, to the variable amount of exposed soil and soil moisture content, to the rapid temporal dynamic 
of the canopy in response to climatic drivers like temperature and water status (Gao et al., 2005) as well 
as to management practices (e.g. planting, harvesting, grazing, crop type, etc.) (Tittebrand et al., 2009). 
Given the large spatial variability of GRA and CRO albedo, future intercomparison of satellite and in 
situ observation for these PFTs should be performed preferentially at selected sites with homogenous 
crop cover and management practices at the MODIS spatial resolution. 

ENF and EBF snow-free albedo do not show any significant seasonal trend and monthly values are 
around 0.1. EBF do not show seasonal variation also at the global scale, while larger variability is 
observed for ENF in winter months. A clear seasonality is shown by DBF+MF and CRO, with higher 
values at the peak of the growing season, in line with the results of other in situ studies (Sellers et al., 
1997). The seasonal pattern of GRA is the opposite, with a winter maximum in albedo. The increase of 
GRA albedo from August, observed both in MODIS and FLUXNET data, is probably due to vegetation 
senescence, with the increase of exposed soil and dead biomass and the drying of stalks and seeds 
during summer months. 
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4. Conclusions 


Surface albedo is a key parameter in the Earth’s energy balance since it affects the amount of solar 
radiation directly absorbed at the planet surface. Its variability in time and space can be retrieved 
through the use of remote sensing products, available nowadays at high temporal and spatial resolution 
from different satellite platforms (e.g. Terra and Aqua MODIS observations are used to retrieve albedo 
every 16 days at a spatial resolution of 500 x 500 m starting March 2000). 

Careful intercomparisons with in situ measurements of surface albedo are crucial to evaluate and 
improve the quality of remote sensing products. In this context the “La Thuile” FLUXNET dataset 
offers an unprecedented opportunity to select sites according to strict conditions of landscape 
homogeneity and across a wide range of PFTs and geographical areas. When compared to the 
latitudinal distribution of surface albedo, the FLUXNET data set clearly shows a bias towards darker 
vegetated areas in the Northern hemisphere. In this respect, the biased distribution of the terrestrial 
dataset in the global albedo domain should be taken into account when using in situ FLUXNET albedo 
data in the parameterization of land surface models. 

The key issue in the intercomparison of satellite and in situ measurements is the spatial 
representativeness of the latter in the satellite pixel. The MODIS pixel can be 100-1000 times larger 
than the footprint of in situ observations, depending on the height of the albedometer above the canopy 
top. For this reason a careful classification of FLUXNET sites has been performed to select only those 
sites with the lower level of heterogeneity. Albedo measurements at the sites passing the selection are 
in excellent agreement for most forest sites, while satellite retrievals underestimate in situ 
measurements for PFTs with larger values of albedo (typically non-forest ecosystems as CRO and 
GRA). 

When seasonal patterns of MODIS albedo are compared, the match is extremely good for all forest 
PFTs. These results quantitatively document the quality achieved with satellite retrievals of surface 
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albedo over a range of forest ecosystems. On the contrary, in non-forest sites satellite retrievals 
systematically underestimate in situ measurements across the whole seasonal cycle with a minimum 
error at the peak of the growing season. However, in the interpretation of these results it should be 
considered that the footprint of the surface radiometric measurements is extremely limited in these 
ecosystems, since sensors are generally installed few meters above the ground. In addition the temporal 
and fine-scale spatial variability of the canopy reflectance can be extremely large in non-forest PFTs 
due to management practices and rapid LAI dynamics (Tittebrand et al., 2009). 

The intercomparison of surface measurements and satellite retrievals is affected by the surface 
heterogeneity at various spatial scale. The sub-pixel variability affects the difference between the signal 
in the footprint of the albedometer and the 500 m MODIS pixel, while the between pixels variability 
affects the accuracy of the inversion algorithm. At smaller scale, the combination of quantitative 
indexes based on variogram analysis of high resolution images confirms that FLUXNET sites in 
herbaceous ecosystems (CRO and GRA) are the less homogenous and should therefore be carefully 
evaluated in point to pixel intercomparisons. At a coarser scale, the variability of the MODIS albedo in 
the 7x7 Km 2 area is related to the mismatch of MODIS retrievals and in situ measurements, further 
corroborating the relevance of large scale landscape homogeneity in the comparison of surface 
measurements and satellite retrievals. 

The results of this investigation clearly show the need for the future to characterize the spatial 
heterogeneity of reference sites using a combination of surface measurements, airborne and finer scale 
satellite imagery (Roman et al., 2009). To overcome the limitations of the point to pixel albedo inter- 
comparisons, future field and airborne surveys should be planned to address the issue of spatial scales 
and landscape heterogeneity. Current approaches for measuring in-situ albedo are not adequate to 
describe mixed or highly heterogeneous landscapes such as mixed forests, open shrublands, savannas 
and croplands. For this purpose only spatially-distributed measurements (i.e., airborne laser scanning, 
network of tall towers) could produce observations at the required scale. For albedo in particular, multi- 
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angle airborne instruments such as AirMISR (Diner et al., 1998), AirPOLDER (Chen et al., 1997), 
NASA's Cloud Absorption Radiometer (Gatebe et al., 2003) and the Compact Airborne Spectrographic 
Imager (CASI) data (Chen et al, 1999) could be effectively used to retrieve surface-level bidirectional 
reflectance measurements at landscape scale and therefore effectively complement moderate resolution 
satellite systems like MODIS and MISR. 
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Figure captions 


Figure 1 - Spatial distribution of the 120 FLUXNET sites for which albedo measurements are 
available. Green dots represent sites selected for the analysis according to plant cover 
homogeneity at 1 Km 2 scale (n=53, visual classification based on high resolution Google 
Earth™ images, e.g. Fig 2). 

Figure 2 - Frequency distribution of synchronous MODIS retrievals and in situ measurements 

classified according to year (panel a) and month (panel b) of observation. The vegetation is 
coded according to the IGBP classification: CRO, croplands; CSF1, closed shrublands; DBF, 
deciduous broad-leaf forests; EBF, evergreen broad-leaf forests; ENF, evergreen needle-leaf 
forests; GRA, grassland; MF, mixed forests; OSF1, open shrublands; SAV, savannas; WET, 
permanent wetlands; WSA, woody savannas. Validation sites are separated into forest (ENF, 
EBF, DBF, MF, SAV and WSA) and non forest (OSH, CSH, CRO, GRA, WET). 

Figure 3 - Global distribution of MODIS albedo (0.05° res.) categorized for latitudinal bands (panel a) 
and Plant Functional Types (panel b). In both panels black dots represent global MODIS 
averages (grey bars show the 1 0-90 percentile intervals) and white dots represent MODIS 
retrievals at the 53 FLUXNET sites reported in Table 1. 

Figure 4 - Examples of high resolution Google Earth™ images used to visually classify the FLUXNET 

2 

sites according to land cover homogeneity. Images cover an area of 1km" area centered at the 
tower coordinates. 

Figure 5 - Time series of synchronous in situ observations (open dots) and MODIS retrievals (black 
dots) at the four sites shown in Fig. 4. 

Figure 6 - MODIS retrievals versus in situ observations of clear-sky and snow-free albedo, a) Averages 
and standard deviations of retrievals (temporal and spatial variability are accounted for) 
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grouped by Plant Functional Types (PFT), b) Average ground measurements and retrievals at 
individual sites classified by PFT. 

Figure 7 - Mean error and mean absolute error of MODIS albedo retrievals at FLUXNET sites. Site 
numbers are reported in Tablel; dot size is proportional to the R 2 of the regression of daily 
MODIS retrievals versus in situ measurements. 

Figure 8 - Mean absolute percentage error of MODIS retrievals versus in situ measurements as a 

function of the standard deviation of MODIS retrievals in the 7x7 Km grid centered at the site. 
Sites located in a more homogeneous landscape (lower spatial variability of albedo) are 
characterized by a lower discrepancy between in situ observations and satellite retrievals. 

Figure 9 - Seasonal average values of clear-sky snow-free albedo at the 53 homogeneous sites as 
measured at the FLUXNET sites or retrieved from MODIS. Winter values (DJF) show the 
lowest correlation. The first panel includes the correlation matrix of the residuals (MODIS 
retrievals- in situ albedo) in the four seasons. 

Figure 10 - Seasonal trends of albedo for the different plant functional types as derived from in situ 
observations (white dots) and MODIS retrievals (black dots). The grey bands show the 
observed variability (10-90 percentiles) of the MODIS retrieved albedo at global scale for each 
plant functional type. 


23 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 


References 


Alton, P. (2009). A simple retrieval of ground albedo and vegetation absoiptance from MODIS satellite 
data for parameterisation of global Land- Surface Models. Agricultural and Forest Meteorology, 149, 
1769-1775. 

Anderson, R.G., Canadell, J.G., Randerson, J.T., Jackson, R.B., Hungate, B.A., Baldocchi, D.D., Ban- 
Weiss, G.A., Bonan, G.B., Caldeira, K., Cao, L., Diffenbaugh, N.S., Gumey, K.R., Kueppers, L.M., 
Law, B.E., Luyssaert, S., & O'Halloran, T.L. (2010). Biophysical considerations in forestry for climate 
protection. Frontiers in Ecology and the Environment, doi: 10.1 890/090 179. 

Baldocchi, D., Falge, E., Gu, L., Olson, R., Hollinger, D., Running, S., Anthoni, P., Bemhofer, C., 
Davis, K., Evans, R., Fuentes, J., Goldstein, A., Katul, G., Law, B., Lee, X., Malhi, Y., Meyers, T., 
Munger, W., Oechel, W., Paw, K.T., Pilegaard, K., Schmid, H.P., Valentini, R., Verma, S., Vesala, T., 
Wilson, K., & Wofsy, S. (2001). FLUXNET: A New Tool to Study the Temporal and Spatial 
Variability of Ecosystem-Scale Carbon Dioxide, Water Vapor, and Energy Flux Densities. Bulletin of 
the American Meteorological Society, 82, 2415-2434. doi: 10.1 175/1520-0477. 

Betts, R.A. (2000). Offset of the potential carbon sink from boreal forestation by decreases in surface 
albedo. Nature, 408, 187-190. 

Betts, R.A., Falloon, P.D., Goldewijk, K.K., & Ramankutty, N. (2007). Biogeophysical effects of land 
use on climate: Model simulations of radiative forcing and large-scale temperature change. 

Agricultural and Forest Meteorology, 142, 216-233. 

Bird, D.N., Kunda, M., Mayer, A., Schlamadinger, B., Canella, L., & Johnston, M. (2008). 
Incorporating changes in albedo in estimating the climate mitigation benefits of land use change 
projects. Biogeosciences Discuss., 5, 1511-1543. 

24 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 


Bonan, G.B. (2008). Forests and Climate Change: Forcings, Feedbacks, and the Climate Benefits of 
Forests. Science, 320, 1444-1449. 

Cescatti, A. (1998). Effects of needle clumping in shoots and crowns on the radiative regime of a 
Norway spruce canopy. Annals of Forest Science, 55, 89-102. 

Chapin, F.S., Randerson, J.T., McGuire, A.D., Foley, J.A., & Field, C.B. (2008). Changing feedbacks 
in the climate-biosphere system. Frontiers in Ecology and the Environment, 6, 313-320. 
doi:10. 1890/080005. 

Chapin, F.S., Sturm, M., Serreze, M.C., McFadden, J.P., Key, J.R., Lloyd, A.FL, McGuire, A.D., Rupp, 
T.S., Lynch, A.FL, Schimel, J.P., Beringer, J., Chapman, W.L., Epstein, FLE., Euskirchen, E.S., 
Hinzman, L.D., Jia, G., Ping, C.L., Tape, K.D., Thompson, C.D.C., Walker, D.A., & Welker, J.M. 
(2005). Role of land-surface changes in arctic summer warming. Science, 310, 657-660. 

Charlson, R.J., Ackerman, A.S., Bender, F.A.M., Anderson, T.L., & Liu, Z. (2007). On the climate 
forcing consequences of the albedo continuum between cloudy and clear air. Tellus Series B-Chemical 
and Physical Meteorology, 59, 715-727. doi: 10.1 1 1 1/j.l 600-0889.2007.00297. 

Chamey, J., Quirk, W.J., Chow, S.FL, & Komfield, J. (1977). A comparative study of the effects of 
albedo change on drought in semi-arid regions. Journal of the Atmospheric Sciences, 34, 1366-1385. 

Chen, J.M., Leblanc, S.G., Cihlar, J.C., Bicheron, P., Leroy, M., Deering, D., & Eck, T. (1997). Studies 
of BRDF in conifer and deciduous boreal forests using the 4-scale model and airborne POLDER and 
ground-based PARABOLA measurements. In, Geoscience and Remote Sensing, 1997. IGARSS '97. 
Remote Sensing - A Scientific Vision for Sustainable Development., 1997 IEEE International (pp. 165- 
167 vol.161) 


25 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 


Chen, J.M., Leblanc, S.G., Miller, J.R., Freemantle, J., Loechel, S.E., Walthall, C.L., Innanen, K.A., & 
White, H.P. (1999). Compact Airborne Spectrographic Imager (CASI) used for mapping biophysical 
parameters of boreal forests. J. Geophys. Res., 104, 27945-27958. 

Chen, Y.M., Liang, S., Wang, J., Kim, H.Y., & Martonchik, J.V. (2008). Validation of MISR land 
surface broadband albedo. International Journal of Remote Sensing, 29, 6971-6983. 
10.1080/01431160802199876. 

Davidson, A., & Wang, S.S. (2004). The effects of sampling resolution on the surface albedos of 
dominant land cover types in the North American boreal region. Remote Sensing of Environment, 93, 
211-224. 10.1016/j.rse.2004.07.005. 

Dickinson, R.E. (1983). Land surface processes and climate surface albedos and energy -balance. 
Advances in Geophysics, 25, 305-353. 

Diner, D.J., Barge, L.M., Bruegge, C.J., Chrien, T.G., Conel, J.E., Eastwood, M.L., Garcia, J.D., 
Hernandez, M.A., Kurzweil, C.G., Ledeboer, W.C., Pignatano, N.D., Sarture, C.M., & Smith, B.G. 
(1998). The Airborne Multi-angle Imaging SpectroRadiometer (AirMISR): instrument description and 
first results. Geoscience and Remote Sensing, IEEE Transactions on, 36, 1339-1349. 

Dirmeyer, P.A., & Shukla, J. (1994). Albedo as a modulator of climate response to tropical 
deforestation. Journal Geophysical Research, 99, 20863-20877. 

Frida A-M, B., Henning, R., Robert J, C., Annica M. L, E., & Norman, L. (2006). 22 views of the 
global albedo - comparison between 20 GCMs and two satellites. Tellus A, 58, 320-330. 

Gao, F., Schaaf, C.B., Strahler, A.H., Roesch, A., Lucht, W., & Dickinson, R. (2005). MODIS 
bidirectional reflectance distribution function and albedo Climate Modeling Grid products and the 


26 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 


variability of albedo for major global vegetation types. Journal of Geophysical Research-Atmospheres, 
110, D01 104. doi: 10.1 029/2004jd005 1 90. 

Gatebe, C.K., King, M.D., Platnick, S., Arnold, G.T., Vermote, E.F., & Schmid, B. (2003). Airborne 
spectral measurements of surface&#821 l;atmosphere anisotropy for several surfaces and ecosystems 
over southern Africa. J. Geophys. Res., 108, 8489. 

Govaerts, Y., & Lattanzio, A. (2008). Estimation of surface albedo increase during the eighties Sahel 
drought from Meteosat observations. Global and Planetary Change, 64, 139-145. 

Hall, A., & Qu, X. (2006). Using the current seasonal cycle to constrain snow albedo feedback in future 
climate change. Geophysical Research Letters, 33. doi:10.1029/2005gl025127. 

Henderson-Sellers, A., & Wilson, M.F. (1983). Surface albedo data for climatic modeling. Review’s of 
Geophysics, 21, 1743-1778. 

Holben, B.N., Eck, T.F., Slutsker, I., Tanre, D., Buis, J.P., Setzer, A., Vermote, E., Reagan, J.A., 
Kaufman, Y.J., Nakajima, T., Lavenu, F., Jankowiak, I., & Smirnov, A. (1998). AERONET— A 
Federated Instrument Network and Data Archive for Aerosol Characterization. Remote Sensing of 
Environment, 66, 1-16. 

Hollinger, D.Y., Ollinger, S.V., Richardson, A.D., Meyers, T.P., Dail, D.B., Martin, M.E., Scott, N.A., 
Arkebauer, T.J., Baldocchi, D.D., Clark, K.L., Curtis, P.S., Davis, K.J., Desai, A.R., Dragoni, D., 
Goulden, M.L., Gu, L., Katul, G.G., Pallardy, S.G., PawU, K.T., Schmid, H.P., Stoy, P.C., Suyker, 
A.E., & Verma, S.B. (2009). Albedo estimates for land surface models and support for a new paradigm 
based on foliage nitrogen concentration. Global Change Biology, 16, 696-710. 


27 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 


Jin, Y., Schaaf, C.B., Woodcock, C.E., Gao, F., Li, X., Strahler, A.H., Lucht, W., & Liang, S. (2003a). 
Consistency of MODIS surface bidirectional reflectance distribution function and albedo retrievals: 2. 
Validation. Journal Geophysical Research, 108. doi:10.1029/2002JD002804. 

Jin, Y.F., Schaaf, C.B., Gao, F., Li, X.W., Strahler, A.H., Lucht, W., & Liang, S.L. (2003b). 
Consistency of MODIS surface bidirectional reflectance distribution function and albedo retrievals: 1. 
Algorithm performance. Journal of Geophysical Research-Atmospheres, 108. 
doi:10.1029/2002jd002803. 

Lawrence, D.M., & Slingo, J.M. (2004). An annual cycle of vegetation in a GCM. Part II: global 
impacts on climate and hydrology. Climate Dynamics, 22, 107-122. doi:10.1007/s00382-003-0367-8. 

Liang, S.L., Fang, H.L., Chen, M.Z., Shuey, C.J., Walthall, C., Daughtry, C., Morisette, J., Schaaf, C., 
& Strahler, A. (2002). Validating MODIS land surface reflectance and albedo products: methods and 
preliminary results. Remote Sensing of Environment, 83, 149-162. 

Liu, J., Schaaf, C., Strahler, A., Jiao, Z., Shuai, Y., Zhang, Q., Roman, M., Augustine, J.A., & Dutton, 
E.G. (2009). Validation of Moderate Resolution Imaging Spectroradiometer (MODIS) albedo retrieval 
algorithm: Dependence of albedo on solar zenith angle. Journal Geophysical Research, 114, D01 106. 

Lucht, W., Hyman, A.H., Strahler, A.H., Barnsley, M.J., Hobson, P., & Muller, J.P. (2000a). A 
Comparison of Satellite-Derived Spectral Albedos to Ground-Based Broadband Albedo Measurements 
Modeled to Satellite Spatial Scale for a Semidesert Landscape. Remote Sensing of Environment, 74, 85- 
98. 

Lucht, W., Schaaf, C.B., & Strahler, A.H. (2000b). An algorithm for the retrieval of albedo from space 
using semiempirical BRDF models. IEEE transactions on geoscience and remote sensing, 38, 977-998. 


28 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 


Masuoka, E., Tilmes, C., Gang, Y., & Devine, N. (2000). Producing global science products for the 
MODerate resolution Imaging Spectroradiometer (MODIS) in MODAPS. In, Geoscience and Remote 
Sensing Symposium, 2000. Proceedings. IGARSS 2000. IEEE 2000 International (pp. 2050-2052 
vol.2055) 

Masuoka, E., Wolfe, R., Sinno, S., Gang, Y., & Teague, M. (2007). A disk-based system for producing 
and distributing science products from MODIS. In, Geoscience and Remote Sensing Symposium, 2007. 
IGARSS 2007. IEEE International (pp. 3043-3046) 

Moody, Eric, G., King, Michael, D., Schaaf, Crystal, B., Platnick, & Steven, G. (2008). MODIS- 
Derived Spatially Complete Surface Albedo Products : Spatial and Temporal Pixel Distribution and 
Zonal Averages. Journal of Applied Meteorology and Climatology, 47, 2879-2894. 

Myhre, G., Kvalevag, M.M., & Schaaf, C.B. (2005). Radiative forcing due to anthropogenic vegetation 
change based on MODIS surface albedo data. Geophysical Research Letters, 32, L21410. 

Ollinger, S.V., Richardson, A.D., Martin, M.E., Hollinger, D.Y., Frolking, S.E., Reich, P.B., Plourde, 

L. C., Katul, G.G., Munger, J.W., Oren, R., Smith, M.L., Paw U, K.T., Bolsta, P.V., Cook, B.D., Day, 

M. C., Martin, T.A., Monson, R.K., & Schmid, H.P. (2008). Canopy nitrogen, carbon assimilation, and 
albedo in temperate and boreal forests: Functional relations and potential climate feedbacks. 
Proceedings of the National Academy of Sciences of the United States of America, 105, 19336-19341. 

Pinty, B., Andredakis, I., Clerici, M., Kaminski, T., Tabemer, M., Verstraete, M.M., Gobron, N., 
Plummer, S., & Widlowski, J.L. (2011a). Exploiting the MODIS albedos with the Two-stream 
Inversion Package (JRC-TIP): 1. Effective leaf area index, vegetation, and soil properties. J. Geophys. 
Res., 116, D09105. 


29 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 


Pinty, B., Tabemer, M., Haemmerle, V.R., Paradise, S.R., Vermote, E., Verstraete, M.M., Gobron, N., 
& Widlowski, J.-L. (2011b). Global-Scale Comparison of MISR and MODIS Land Surface Albedos. 
Journal of Climate, 24, 732-749. doi: 1 0.1 175/20 10JCLI3709.1. 

Pirazzini, R. (2004). Surface albedo measurements over Antarctic sites in summer. J. Geophys. Res., 
109, D20118. 

Pirazzini, R., Vihma, T., Granskog, M.A., & Cheng, B. (2006). Surface albedo measurements over sea 
ice in the Baltic Sea during the spring snowmelt period. Annals of Glaciology, 44, 7-14. 

Randerson, J.T., Liu, H., Flanner, M.G., Chambers, S.D., Jin, Y., Hess, P.G., Pfister, G., Mack, M.C., 
Treseder, K.K., Welp, L.R., Chapin, F.S., Harden, J.W., Goulden, M.L., Lyons, E., Neff, J.C., Schuur, 
E.A.G., & Zender, C.S. (2006). The Impact of Boreal Forest Fire on Climate Warming. Science, 314, 
1130-1132. doi: 10.1 126/science. 1 132075. 

Roman, M.O., Schaaf, C.B., Lewis, P., Gao, F., Anderson, G.P., Privette, J.L., Strahler, A.H., 
Woodcock, C.E., & Barnsley, M. (2010). Assessing the coupling between surface albedo derived from 
MODIS and the fraction of diffuse skylight over spatially-characterized landscapes. Remote Sensing of 
Environment, 114, 738-760. 

Roman, M.O., Schaaf, C.B., Woodcock, C.E., Strahler, A.H., Yang, X., Braswell, R.H., Curtis, P.S., 
Davis, K.J., Dragoni, D., Goulden, M.L., Gu, L., Hollinger, D.Y., Kolb, T.E., Meyers, T.P., Munger, 
J.W., Privette, J.L., Richardson, A.D., Wilson, T.B., & Wofsy, S.C. (2009). The MODIS (Collection 
V005) BRDF/albedo product: Assessment of spatial representativeness over forested landscapes. 
Remote Sensing of Environment, 113, 2476-2498. 

Rotenberg, E., & Yakir, D. (2010). Contribution of Semi-Arid Forests to the Climate System. Science, 
327, 451-454. doi: 10. 1 126/science. 1179998. 


30 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 


Running, S.W., Loveland, T.R., Pierce, L.L., Nemani, R.R., & Hunt, E.R. (1995). A remote sensing 
based vegetation classification logic for global land cover analysis. Remote Sensing of Environment, 

51, 39-48. 

Salomon, J., G., Schaaf, C., B., Strahler, A., H., Feng, G.A.O., & Yufang Jin, G. (2006). Validation of 
the MODIS bidirectional reflectance distribution function and albedo retrievals using combined 
observations from the Aqua and Terra platforms. IEEE transactions on geoscience and remote sensing, 
44, 1908-1925. 

Salomonson, V.V., Barnes, W.L., Maymon, P.W., Montgomery, H.E., & Ostrow, H. (1989). MODIS - 
advanced facility instrument for studies of the earth as a system. IEEE transactions on geoscience and 
remote sensing, 27, 145-153. 

Schaaf, C.B., Gao, F., Strahler, A.H., Lucht, W., Li, X.W., Tsang, T., Strugnell, N.C., Zhang, X.Y., Jin, 
Y.F., Muller, J.P., Lewis, P., Barnsley, M., Hobson, P., Disney, M., Roberts, G., Dunderdale, M., Doll, 
C., d'Entremont, R.P., Hu, B.X., Liang, S.L., Privette, J.L., & Roy, D. (2002). First operational BRDF, 
albedo nadir reflectance products from MODIS. Remote Sensing of Environment, 83, 135-148. 

Sellers, P.J., Hall, F.G., Kelly, R.D., Black, A., Baldocchi, D., Berry, J., Ryan, M., Ranson, K.J., Crill, 
P.M., Lettenmaier, D.P., Margolis, H., Cihlar, J., Newcomer, J., Fitzjarrald, D., Jarvis, P.G., Gower, 
S.T., Halliwell, D., Williams, D., Goodison, B., Wickland, D.E., & Guertin, F.E. (1997). BOREAS in 
1997: Experiment overview, scientific results, and future directions. Journal of Geophysical Research- 
Atmospheres, 102, 28731-28769. 

Stroeve, J., Box, J.E., Gao, F., Liang, S., Nolin, A., & Schaaf, C. (2005). Accuracy assessment of the 
MODIS 16-day albedo product for snow: comparisons with Greenland in situ measurements. Remote 
Sensing of Environment, 94, 46-60. 


31 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 


Susaki, J., Yasuoka, Y., Kajiwara, K., Honda, Y., & Hara, K. (2007). Validation of MODIS Albedo 
Products of Paddy Fields in Japan. Geoscience and Remote Sensing, IEEE Transactions on, 45, 206- 
217. 

Tian, Y., Dickinson, R.E., Zhou, L., Myneni, R.B., Friedl, M., Schaaf, C.B., Carroll, M., & Gao, F. 
(2004). Land boundary conditions from MODIS data and consequences for the albedo of a climate 
model. Geophysical Research Letters, 31, L05504. 

Tittebrand, A., Spank, U., & Bemhofer, C.H. (2009). Comparison of satellite- and ground-based NDVI 
above different land-use types. Theoretical and Applied Climatology, 98, 171-186. 

Vavrus, S., Ruddiman, W.F., & Kutzbach, J.E. (2008). Climate model tests of the anthropogenic 
influence on greenhouse-induced climate change: the role of early human agriculture, industrialization, 
and vegetation feedbacks. Quaternary’ Science Reviews, 27, 1410-1425. 

Wang, K., Liang, S., Schaaf, C.L., & Strahler, A.H. (2010). Evaluation of Moderate Resolution 
Imaging Spectroradiometer land surface visible and shortwave albedo products at FLUXNET sites. 
Journal Geophysical Research, 115, D17107. 


32 


Table 1 - Site characteristics derived from the FLUXNET database. N obs indicates the number of days of 
synchronous recordings of MODIS and in situ albedo. The averages of MODIS retrievals and in situ 
observations are reported in the two rightmost columns. 


N 

SITE ID 

COUNTRY 

PFT 

Lat. 

[deg] 

Long. 

[deg] 

N Obs 

MODIS 

albedo 

In situ 
albedo 

ST score 

1 

AU Turn 

Australia 

EBF 

-35.66 

148.15 

733 

0.11 

0.11 

9.75 

2 

AU Wac 

Australia 

EBF 

-37.43 

145.19 

275 

0.09 

0.10 

11.30 

3 

BR Cax 

Brazil 

EBF 

-1.72 

-51.46 

67 

0.12 

0.12 

5.33 

4 

BR Sa3 

Brazil 

EBF 

-3.02 

-54.97 

79 

0.13 

0.12 

6.23 

5 

BW Ghg 

Botswana 

SAV 

-21.51 

21.74 

28 

0.16 

0.18 

0.67 

6 

BW Ghm 

Botswana 

WSA 

-21.20 

21.75 

28 

0.18 

0.17 

NA 

7 

BW Mai 

Botswana 

WSA 

-19.92 

23.56 

252 

0.16 

0.14 

5.42 

8 

CA Cal 

Canada 

ENF 

49.87 

-125.33 

580 

0.09 

0.09 

5.99 

9 

CA Ca3 

Canada 

ENF 

49.53 

-124.90 

552 

0.13 

0.14 

0.64 

10 

CA NS6 

Canada 

OSH 

55.92 

-98.96 

353 

0.10 

0.12 

NA 

11 

CASF2 

Canada 

ENF 

54.25 

-105.88 

363 

0.11 

0.11 

1.72 

12 

CASF3 

Canada 

ENF 

54.09 

-106.01 

362 

0.10 

0.10 

3.14 

13 

CAWP1 

Canada 

MF 

54.95 

-112.47 

353 

0.11 

0.13 

1.65 

14 

CZ BK1 

Czech Republic 

ENF 

49.50 

18.54 

148 

0.09 

0.10 

3.18 

15 

DE Geb 

Germany 

CRO 

51.10 

10.91 

328 

0.17 

0.18 

0.90 

16 

DE Hai 

Germany 

DBF 

51.08 

10.45 

451 

0.13 

0.13 

1.97 

17 

DE Kli 

Germany 

CRO 

50.89 

13.52 

256 

0.16 

0.19 

0.62 

18 

DE Tha 

Germany 

ENF 

50.96 

13.57 

477 

0.10 

0.07 

5.99 

19 

DE Wet 

Germany 

ENF 

50.45 

11.46 

375 

0.07 

0.05 

1.55 

20 

ES ES2 

Spain 

CRO 

39.28 

-0.32 

298 

0.14 

0.13 

0.91 

21 

ES LMa 

Spain 

SAV 

39.94 

-5.77 

455 

0.16 

0.19 

1.89 

22 

FR Fon 

France 

DBF 

48.48 

2.78 

162 

0.14 

0.13 

0.51 

23 

FR Hes 

France 

DBF 

48.67 

7.07 

474 

0.15 

0.14 

1.44 

24 

FR Pue 

France 

EBF 

43.74 

3.60 

401 

0.11 

0.12 

0.87 

25 

GF Guy 

French Guyana 

EBF 

5.28 

-52.93 

216 

0.12 

0.10 

1.92 

26 

HU Bug 

Hungary 

GRA 

46.69 

19.60 

523 

0.16 

0.20 

1.74 

27 

IE Dri 

Ireland 

GRA 

51.99 

-8.75 

160 

0.20 

0.22 

1.05 

28 

IT Bon 

Italy 

ENF 

39.48 

16.54 

174 

0.10 

0.08 

2.43 

29 

IT Col 

Italy 

DBF 

41.85 

13.59 

183 

0.14 

0.15 

2.60 

30 

ITSRo 

Italy 

ENF 

43.73 

10.28 

512 

0.08 

0.09 

3.10 

31 

JP Mas 

Japan 

CRO 

36.05 

140.03 

177 

0.12 

0.13 

1.15 

32 

KR Kwl 

Korea 

MF 

37.75 

127.16 

216 

0.09 

0.09 

3.81 

33 

NL Cal 

Netherlands 

GRA 

51.97 

4.93 

369 

0.17 

0.22 

1.03 

34 

NL Lan 

Netherlands 

CRO 

51.95 

4.90 

108 

0.17 

0.18 

1.44 

35 

NL Loo 

Netherlands 

ENF 

52.17 

5.74 

404 

0.11 

0.09 

1.40 

36 

PT Esp 

Portugal 

EBF 

38.64 

-8.60 

437 

0.13 

0.11 

0.86 

37 

RU Che 

Russia 

MF 

68.61 

161.34 

104 

0.15 

0.16 

0.54 

38 

SE Nor 

Sweden 

ENF 

60.09 

17.48 

234 

0.10 

0.09 

3.84 

39 

UKGri 

UK 

ENF 

56.61 

-3.80 

15 

0.10 

0.09 

1.48 

40 

US Aud 

USA 

GRA 

31.59 

-110.51 

1244 

0.18 

0.21 

0.62 

41 

US Bnl 

USA 

ENF 

63.92 

-145.38 

104 

0.11 

0.09 

0.56 

42 

US Bol 

USA 

CRO 

40.01 

-88.29 

721 

0.16 

0.19 

NA 

43 

US Bo2 

USA 

CRO 

40.01 

-88.29 

289 

0.16 

0.19 

NA 

44 

US Fmf 

USA 

ENF 

35.14 

-111.73 

247 

0.10 

0.12 

3.08 

45 

US FPe 

USA 

GRA 

48.31 

-105.10 

775 

0.15 

0.16 

1.05 

46 

US Fuf 

USA 

ENF 

35.09 

-111.76 

231 

0.10 

0.13 

1.72 

47 

US IB1 

USA 

CRO 

41.86 

-88.22 

349 

0.15 

0.16 

0.76 

48 

US Ivo 

USA 

WET 

68.49 

-155.75 

62 

0.18 

0.19 

NA 

49 

US MMS 

USA 

DBF 

39.32 

-86.41 

639 

0.13 

0.13 

8.75 

50 

US MOz 

USA 

DBF 

38.74 

-92.20 

563 

0.12 

0.11 

2.17 

51 

US SRM 

USA 

WSA 

31.82 

-110.87 

826 

0.17 

0.16 

2.13 

52 

US WCr 

USA 

DBF 

45.81 

-90.08 

487 

0.14 

0.15 

1.82 

53 

ZA Kru 

South Africa 

SAV 

-25.02 

31.50 

447 

0.15 

0.15 

1.34 


Table 2 - PFT averages of mean albedo (ground observation and satellite retrievals) and geostatistical indexes 
of landscape heterogeneity based on Landsat images. The analysis is limited to the 48 sites with available 
Landsat7 images. Records are shown in descending order of ST score . 


PFT 

N sites 

MODIS 

In-situ 

MAE 

> 

u 

Rse 

Rsv 

Rst 

ST score 

EBF 

7 

0.115 

0.109 

0.014 

17% 

20% 

-3% 

48% 

5.18 

WSA 

2 

0.167 

0.148 

0.020 

14% 

52% 

31% 

-3% 

3.77 

DBF 

7 

0.138 

0.136 

0.013 

42% 

15% 

78% 

10% 

2.75 

ENF 

15 

0.100 

0.095 

0.017 

24% 

23% 

44% 

17% 

2.65 

MF 

3 

0.115 

0.126 

0.016 

41% 

39% 

41% 

56% 

2.00 

SAV 

3 

0.158 

0.174 

0.017 

9% 

58% 

61% 

-13% 

1.30 

GRA 

5 

0.172 

0.203 

0.033 

13% 

82% 

63% 

-1% 

1.10 

CRO 

6 

0.155 

0.164 

0.020 

11% 

85% 

7% 

24% 

0.96 
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